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Abstract 

Recent evidence from observation of the flavor-changing neutral current decay B® — > [i + fJ,~ by the 

Qh 

O , LHCb collaboration B(Bg — > fi + u~) = (3.2jl^'|) x 10~ 9 is consistent with the latest standard model 
D 

prediction Bsm(B^ s ~~ ^ / i+ A i ) = (3.23 ± 0.27) x 10 9 . While new physics can still affect this decay 

amplitude, its contribution is certainly not the dominant one. We analyze branching ratios of the decays 
\q B® — > u + [i~X, with X = 7 or vv, which can mimic B® — > u + u~ on portions of the parameter space where 
t-h X is soft. We perform a model- independent standard model calculation of those processes incorporating 

heavy quark and chiral symmetries of QCD. We show that the considered contributions contribute to 

Bg — > u + fi~ at a sub-percent level. 
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I. INTRODUCTION 



The rare leptonic decay of the B® into a dimuon pair, B® — > is an example of a flavor- 

changing neutral current (FCNC) process. Studies of such decay processes not only play an 
important role in determining electroweak and strong interaction parameters of the standard model 
(SM) of particle physics, but also serve as sensitive probes of possible physics beyond the standard 
model While recent evidence for observation of B° s — > from LHC-b collaboration [2] 

, as well as an earlier result from CDF [3] preclude any spectacular new physics (NP) effect, 
there is still room for NP to influence this decay. It is then important to have a firm evaluation of 
B(B® — > in the SM [HE] and a firm understanding that experimentally-observed branching 

ratio 

&LBCb{B° s -+ /xV) = (3.2j£|) x 1(T 9 

b C df(b° s -> ^V) = (i-8±Si) x ict 8 (1) 

actually corresponds to the B° s — > /i + /i~ transition. 

It is well known that the B° s — > ^ "jjT decay is helicity suppressed in the SM by rnf l /m 2 Bg due 
to the left handed nature of weak interactions [6]. This effect arises from the necessary spin flip 
on the outgoing back-to-back lepton pair in order to conserve angular momentum since the initial 
state meson is spinless. 

This suppression is absent in B° s decays where the muon pair is produced with one or more 
additional particles in the final state that can carry away a unit of angular momentum, such as 
B Q S — > ^ or Bg — > ^fi^u^u^. This means that, in general, those processes could have sizable 
total branching ratios, comparable to that of B° s — > despite being suppressed by other 

small parameters (such as ocem for B° s — > /i + /i~7) |7]. If, in addition, the final state photon or 
vv is undetected, while the invariant mass of fi + [i~ pair is close to m#o, then the experimentally- 
measured branching ratio would correspond to 

B exp (B° s -> //+/0 = B(B° S -> 



B(B° -> p+irX) 



\m(ij.+ fj,-)&m Bs ) 



X 



B(B° M + A i-) 



(2) 
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where X is an undetected particle or a group of particles. The contribution of B® — > jj, + /2^X 
would depend on how well X could be detected in a particular experiment, as well as on whether 
B° s — > fi + fi~X has any kind of resonance enhancement that is not well modeled by background 
models chosen by a particular experiment in a given window of m(/i + /i~), as well as the size of 
that window. For example, for X = 7, most current searches use di-lepton energy cuts that would 
correspond to an allowable soft photon of up to 60 MeV. For B — > tvn transition and X = 7 
similar effects were discussed in [8HTU]. and for X being light particles - in [11]. In the following 
we shall concentrate on the amplitudes that are non- vanishing in the m M — > limit. 

II. B° /i+/x~7 TRANSITION 

Due to higher backgrounds in hadron collider experiments soft photons emitted in B® — > £ + £~'~f 
could be hard to detect, so this background could be quite important. This decays were previously 
analyzed in [7J , where a form-factor model-dependent calculation was performed [12] . The analysis 
presented in [7] was mainly geared towards kinematical regimes where the photon is sufficiently 
hard to be detected; in fact, low-energy cut-offs were introduced on photon energies. We apply a 
model-independent approach that incorporates both heavy quark symmetry for hadrons containing 
a heavy quark with mass rriQ » Aq C d, and chiral SU(3)l x SU(3)r symmetries in the m q — > 
limit [T31 [Ti] . We organize our calculations in terms of an expansion in 1/m^ and examine the 
contribution of terms up to leading order in 0(l/rrib)- 

Similarly to B — > Ivf) [9], the decay amplitude for B s — > / u + / u _ 7 transition can be broken into 
two generic parts containing internal bremsstrahlung (IB) and structure dependent (SD) contri- 
butions. The bremsstrahlung contributions are still helicity suppressed, while the SD contribution 
contain the electromagnetic coupling a but are not suppressed by the lepton mass. Phenomeno- 
logically, the origin of that can be understood as follows. When the soft photon in emitted from 
the B s meson, heavy intermediate states including the J p = 1~ B* vector meson state become 
possible. This lifts helicity suppression since the lepton pair couples directly to the spin 1 meson. 
In the kinematic regime where the photon is soft, we expect that the significant contribution comes 
only from the vector B* resonance for reasons analogous to the B* pole dominance in B — > -niv at 
near zero pionic recoil energies [15J. This is because in the large m;, limit the B* and B s become 



3 




FIG. 1: Resonant contributions to -B, 



degenerate and the residual mass splitting is mg. —vtlb ~ 0(1/ mi,) [S]. Therefore the excitation 
of the B* does not require much energy. There are two diagrams containing an intermediate B* 
as seen in Fig. [TJ In the kinematic region of interest where E 1 < 60MeV, Fig. [T] (b) where the 
intermediate B* decays to an on-shell soft photon is (l/M B o) suppressed and will be neglected. 
Similarly, a contribution of Fig. [I] (d) is formally (l/M#o), so it will be neglected in what follows. 

The calculation of soft photon effects should carefully deal with soft divergencies. Those are 
cancelled between one-loop radiative corrections to B s — > fi + fi~ and B s — > fi + fi~^f. 

We employ Heavy Meson Chiral Perturbation Theory (HMxPT) to calculate Fig. [I] (a). The 
heavy meson superfield H a contains both the B Q S and B*° bosons, 



1+ 



(V-5 a75 ), 



H a = 7 (i/ a ) 7 



UP 



(3) 



where the indices a and b reflect the light quark flavor indices. The light mesons are introduced 
through the matrix = + where the field £ = exp(ill//) is defined in terms of a 

3x3 unitary matrix containing the octet of pseudo-Goldstone bosons 



V2 



n 



= + 




K+ \ 






K° 


K~ 


K° 





(4) 
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FIG. 2: One loop corrections to ji. The double lines denote the heavy mesons B and B* while the single 
line denotes the goldstone bosons 

To evaluate diagram Fig. [I] (a) we need an amplitude for a5-> B*^ transition as 

= M b^ + »- x x M« Bs ^ Btv (5) 

B* 

The amplitude for B — > 5*7 is conventionally parameterized as 

M Bs ^b* s1 = -iefxrjlv^ete^, (6) 

where k is the 4-momentum of the photon, v the velocity of the decaying heavy meson, rj is the 
vector meson polarization, and e is the photon polarization. The strength of the transition is 
described by the magnetic moment, /i, which receives contributions from the photon coupling to 
both heavy and the light quark components of the electromagnetic current [T7] , 

H = fj, b + ii £ . (7) 

The bottom quark contribution is fixed by heavy quark symmetry to be \i b = Q b /m b = — 1/ (3m b ), 
while the light quark contribution can be computed, to one loop, in HM^PT. The relevant effective 
Lagrangian is [T71 HE] 

C p = ^Tr(H a H b a^F, u Qt) + | Tr {H a H bl ^d^ - id^) ba ) , (8) 

where Tr is a trace over the Dirac indices, and /3 is a coupling constant parameterizing a local 
contribution to the light quark magnetic moment. We include the most important one-loop 
correction, which is shown in Fig. [2] 
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The effective magnetic moment for the B s — > 5*7 transition is then 



3m b 3^ ' " 47r/|' 



(9) 



where g is the coupling constant, and ttlk, fx are the mass and decay constant of the Kaon 
respectively. The constants (3 and g can be extracted from a combination of the experimental D* + 
branching ratios, B(D* + ->- D+ 7 ) = 0.016 ± 0.004 and B(D* + ->- D°n + ) = 0.677 ± 0.005, and 
the total width, where the newest preliminary result from BaBar collaboration is reported to be 
T D *+ = 83.5 ± 1.7 ± 1.2 KeV [T9]. The decay widths for these processes using the method above 
are given by 



F(D* + -»■ L> + 7 ) 



7T 



"L .v / _2 1 2 m ^ 

3m c 3^ + ^4tt/2 



(10) 



'ir 



This yeilds the approximate values of the coupling constants, g ~ 0.552 and ~ 7.29GeV 1 . 
With Eqgthis gives us \fi eff \ « 1.13 GeV" 1 . 

To complete evaluation of Fig. [l] (a) in Eq. d5J), we evaluate the 5* — >■ transition. The 

effective hamiltonian describing the weak b — > s£ + £~ transition is 



H 



G, 



b->see 



--Vt b V ts 



v/2 8tt 2 [ 



(12) 



where q u = (pi+ + pi-)„ is the momentum of the lepton pair and Ci are scale-dependent Wilson 
coefficients. The matrix element for B* — > fip, is then 



G 



(13) 
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where 7?^ and v^ 1 are the polarization and 4-velocity of the vector meson respectively. We defined 
(0\s L ^b L \B*) = 77 M / B j/2, and (0|sa^(l + ^)b\B* s ) = M B f Bs [ie^v a r] + v»rf - v v rf\, with 
f B * = MbJbs ED]- This gives for the amplitude of Fig. [T| (a) 



G 



x 



y/2 w ts ^ en q 2 -M 2 Bt 
(2C 7 m b - C 9 M B ») [u Pl -f u v P2 ] - C w M B *[u Pl -f u 'y 5 v P2 



The other contribution that is leading the M Bs — > 00 limit is given in Fig. [T] (c) 



(14) 



M 



Mho^m x x M 
mi 



(15) 



Employing vector-meson dominance, and using the definition of the vector meson decay con- 
stant (O|s7 M s|0) = f<f,m<f>r)%, where 77^ is the polarization of the (ft meson, and (7|s(— ieQ s /A)s\<f>) = 
(-teQ s )e;(0\srs\<P), 

M1^, = l -eUm^el. (16) 



Again, we calculate M.^^,^ using (HM^pT). For the short distance contributions we use the 



effective Hamiltonian describing b — > sit transitions in Eq. 12, as well as the effective Hamiltonian 
for b — >■ sj, 



s~1 

n b -> sl = -^V tb V t * a ^m b C 7 ^ 2 ) -sa^(l + 75)6 ■ F^. 



(17) 



In order to bosonize the quark currents found in Eqs. (12) and (17) we introduce the light vector 
octet to the HM^pT [20], 



■ 9v 
V2 



_|_ JfL 

K*- 



K* 
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The bosonized currents S7 M (1 — 75)6 and sa^il + 75)6 are, respectively, 
Lt=a x ( l5 H b {ff)Ul), 

IZ^i\sT^ ~ \^\( lh H b [ la {p p ) bc - 7/»WfcJ&>. (19) 

A numerical value of cti = — 0.07 GeV 1 ^ 2 [20] will be used for our calculations. Keeping only the 
gauge invariant portion, the amplitude for the decay with an intermediate 0(1020) is 



A{B° S fifKj) -»■ /i/i 7s ) = G F V tb V t 



e 3 hg (j> aiC 7 m b 



e a\ ft ' (Pi + P2))[Upi 



2A^^fM B ~m dp (p 1 -p 2 ) 
(Pi + P2T[u P1 fiv P2 ] + %^ va{i k a {j) X + p 2 )fi\u v ^ v v 



(20) 



We checked that other contributions to the decay are smaller then the ones considered above. We 
considered the bremsstrahlung diagrams where a soft photon is emitted from one of the outgoing 
leptons. These diagrams will result in an infrared divergence in the soft region, which has been 
shown to cancel with the 1-loop QED vertex corrections [21] . The vertex corrections, as well as the 
bremsstrahlung contributions, will remain suppressed by a power of the lepton mass. Therefore 
the remaining non-divergent contributions from both the bremsstrahlung and vertex corrections 
to final states with either an electron or a muon would not be significant. 

The only contribution to the ampltitude from the effective hamiltonian describing the weak 



transition in Eq.(12) ends up being the 0\§ operator. This come from obtaining the matrix 



elements for the pseudoscalar meson, 



<0|(5y*(l - 7S)6)IS a > i/sPg, 

(0|(^(1+7 5 )&)|£Q = 0, 



(21) 
(22) 



where is the decay constant of the B meson. With these definitions and using the conservation 
of the vector current we can arrive at an expression for the amplitude 



M B rem = ie — — =— V tb V t J B C W m ^ 

2y 2 



4h 



B 



Ptt- -k p^+ -k 



75^ 



(23) 
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where e M and k are the polarization and momentum of the photon respectively. Just as we would 
expect from the helicity structure involved, the amplitude for the bremsstrahlung contribution is 
proportional to the lepton mass. So in the limit m^ — > 0, this contribution should be negligible 
compared to the non helicity-suppressed contributions. 

Putting everything together, the distribution of the decay width as a function of the kinematic 
variable s = (P Bs - k) 2 /M B = q 2 /M B , in the limit mi 0, 



dT 

ds 



dT 

ds 



+ 



dT 

ds 



+ 



dT 

ds 



(24) 



where the decay distributions are given for the two different resonance amplitudes and their 
interference. 



dT 

ds 

dT 

ds 

dT 

ds 



B: 



XckmMIjI^j [(\C 9 \ 2 + \C XQ \ 2 )x B% + AC 2 x b -4C 7 C 9 x b x B *] 
lQC 2 f 2 g 2 m 2 a 2 



8(1 



{s-x 2 ) 2 ' 



CRM 



CRM 



9 ml 



(25) 



(C 7 C 9 x B * - 2C 2 7 x b ) 



:i-s) 



s — X 



2 ' 



where we have defined Xqrm = (G 2 F \V tb V t *\ 2 M Bs a EM )/(768ii 4 '), x b = m b /MB s , and xb* = 
M B */M Bs . We use the Wilson coefficients Cj(A) choosing the scale at A ~ m b ~ 5GeV, 
with C 7 = 0.312, C 9 = -4.21 and C w = 4.64 [7] [22]. The CKM matrix elements are 
|V t6 V^| = (4.7 ± 0.8) x 10~ 2 [23J . With the most recent lattice calculation of f Bs is w 228 



MeV [24J. Note that, when integrated over the endpoint window the last two terms in Eq. (24) 
are much smaller than the first one. The interference contribution is destructive and is 



B(B S ii + n 7ij<6o) 
B{B S ->■ /i + / u _ 7 S<30 o) 



-5.0 x 10 



1.1 x 10" 



-17 



-14 



(26) 
(27) 



which are both much smaller than the B* contribution alone. 

The normalized differential spectrum in s is shown in Fig. ([3]). The photon energy is related to 
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FIG. 3: Normalized differential spectrum in s. The grey shaded region corresponds to the contribution 
from a soft photon energy cut at E~ ~ 60 MeV. 

the invariant mass as E-y = (1 — s)Mb/2, so we can integrate the differential spectrum over the 
required corresponding kinematic region in photon energy to obtain the decay width. 



Integrating Eq.(24) over the kinematic region corresponding to a soft photon cut of E~ ~ 60, 300 



MeV we get the respective branching ratios 

B(B S -> /i V7£<6o) = 1-6 x 10~ 12 , (28) 
B{B a -> // + /i-7E<30o) = 1.1 x 10- 10 , (29) 

which are quite too low to affect experimental determination of the branching ratio B s — > fi + fi~ , 
agreeing with the estimates of Ref. @] where B SM (B° S -»■ /J, + fM~) = (3.23 ± 0.27) x 10" 9 . 



III. B° s -> n+n-i/pPp TRANSITION 

Because of the Glashow-Illiopulous-Maiani (GIM) mechanism, the SM loop diagram for the 
helicity-suppressed B® — >■ decay is dominated by the intermediate top quark despite being 

suppressed by the CKM factors |V^V^*| 2 . A transition similar to the ones described above, which on 
a portion of the available phase space looks like B® — > is the tree- level decay B° s — > fi + {i~vv. 

The dominant tree-level contribution for this process is depicted in Fig. |4} This decay can 
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FIG. 4: B s -»• n+yTvv 

have a contribution to the background, which appears only below q 2 = Mg g and, if numerically 
significant, can affect the extraction of B(B S — > ^r). This process is neither loop-dominated 
nor is it helicity suppressed. It nevertheless has a kinematic phase space suppression due to the 
four-particle final state. For the B s meson decay, an intermediate charm quark will give the 
largest contribution since the intermediate top quark diagram will be suppressed by the mass of 
the top quark. Also, the up quark contribution is suppressed by V u t,V* s ~ A 4 whereas the charm 
contribution is only suppressed by V c hV* s ~ A 2 , where A ~ 0.22. 
The transition amplitude for this process is simple, 

Mb^-ub = ^WS<0|*y„(l - T5) # ± ^ 7a(l - j 5 )b\B s )L^ 2 , (30) 
z p c m c 

where L a = /i7 a (l — 75)^. In the rest frame of the decaying meson we can reduce the phase 
space integral's dependence to five independent Lorentz invariants. In the same fashion as in [25] 
we define these invariants as 

5*12 = ijPix- + PtM+f, S 13 = (p M - + p p ) 2 , S 3 4 = (p p + p u ) 2 , 

5*123 = ijPtx- + P^+ + Puf , 5*134 = (Pf,- +p P + p u ) 2 . (31) 



Our width then becomes 



dF = W / (ssO ^W^l 2 dSl2dSmdSl3dSlM , (32) 
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where A4 is the symmetric Gram determinant 



A 4 (q,r,s,t) 



q q ■ r q ■ s q ■ t 
r ■ q r 2 r ■ s r ■ t 
2 s-t 



s ■ q s ■ r s 
t ■ q t ■ r t ■ s 



t 2 



(33) 



In order to avoid the divergence of l/(— A 4 ) 1//2 on the boundary, suitable variable changes can be 
made thereby making the singularity integrable. We define 



S 



134 



6 + sin(S , i 34 )(6 2 -4ac) 



1/2 



2a 

S 13 = 4(-a) 1 / 2 S 13 + m 2 



(34) 



where a, b and c are the parameters solved by 



aS 134 + ^134 + c . 



(35) 



The limits of integration are calculated in [23], resulting in our partial width 



dY 



dS 12 (4vr)6M3 J Si 



M' 2 



/•l pir/2 

dS U3 I dS 3A / dSi3 / dSiu \M Bs ^+ 

Jmj/S 12 J-tt/2 



(36) 



where £ = (M 2 — 6*123) (5*123 — 5*i2)/'S'i23- We define the cut on missing energy as S^ t (E cut ) = 
M 2 — 2M(E cut ) which gives us a lower limit on S12 for the final integral in order to obtain the decay 
width. The branching ratios for this contribution can then calculated using numerical phase-space 
integration for various cuts including the one that corresponds to the invariant mass range seen 
at the LHCb 0. 



BR [B. ->• fi + n~vv] 



E cut =60MeV 



BR \B S ^ iJ+irvV] 



:300MeV 



1.6 x 10" 25 



1.4 x 10~ 18 . 



(37) 



As we can see, the due to enormous phase space suppression (we are only interested in a small sliver 
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of the available four-particle final state), the possible contribution from this decay is unimportant 
for experimental analyses. 

IV. CONCLUSIONS 

Recent evidence of observation of the flavor-changing neutral current decay B® — > fi + ^ by 
the LHCb collaboration B(B° S — > = (3.2^2) x 10~ 9 is consistent with the latest standard 

model predictions. We analyzed branching ratios of the decays B® — > [x + fi~X, with X = 7 or uu, 
which can mimic B° s — > fi + fi~ on portions of the parameter space where X is soft. We performed a 
model-independent calculation of those processes incorporating heavy quark and chiral symmetries 
of QCD. Our calculations concentrated on the contributions that are not helicity suppressed by 
powers of and leads to a correction to the SM prediction of approximately 3% at a photon 
energy cut of 300 MeV and less than 1% at a cut of 60 MeV from soft photon contributions to 
the decay B s — > fi + n~ . The possible contamination from B s — > fi + {i~vv is even smaller, at the 
sub percent level. 

We would like to thank Gil Paz and Rob Harr for useful discussions. This work was supported 
in part by the U.S. Department of Energy under Contract |de-sc 0007983, 
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